Genomic epidemiology offers high resolution estimates of serial intervals for COVID-19.

Serial intervals - the time between symptom onset in infector and infectee - are a fundamental quantity in infectious disease control. However, their estimation requires knowledge of individuals' exposures, typically obtained through resource-intensive contact tracing efforts. We introduce an alternate framework using virus sequences to inform who infected whom and thereby estimate serial intervals. We apply our technique to SARS-CoV-2 sequences from case clusters in the first two COVID-19 waves in Victoria, Australia. We find that our approach offers high resolution, cluster-specific serial interval estimates that are comparable with those obtained from contact data, despite requiring no knowledge of who infected whom and relying on incompletely-sampled data. Compared to a published serial interval, cluster-specific serial intervals can vary estimates of the effective reproduction number by a factor of 2-3. We find that serial interval estimates in settings such as schools and meat processing/packing plants are shorter than those in healthcare facilities.

Proficiency testing for SARS-CoV-2 whole genome sequencing.

Extensive studies and analyses into the molecular features of severe acute respiratory syndrome related coronavirus 2 (SARS-CoV-2) have enhanced the surveillance and investigation of its clusters and transmission worldwide. The whole genome sequencing (WGS) approach is crucial in identifying the source of infection and transmission routes by monitoring the emergence of variants over time and through communities. Varying SARS-CoV-2 genomics capacity and capability levels have been established in public health laboratories across different Australian states and territories. Therefore, laboratories performing SARS-CoV-2 WGS for public health purposes are recommended to participate in an external proficiency testing program (PTP). This study describes the development of a SARS-CoV-2 WGS PTP. The PTP assessed the performance of laboratories while providing valuable insight into the current state of SARS-CoV-2 genomics in public health across Australia. Part 1 of the PTP contained eight simulated SARS-CoV-2 positive and negative specimens to assess laboratories' wet and dry laboratory capacity. Part 2 involved the analysis of a genomic dataset that consisted of a multi-FASTA file of 70 consensus genomes of SARS-CoV-2. Participating laboratories were required to (1) submit raw data for independent bioinformatics analysis, (2) analyse the data with their processes, and (3) answer relevant questions about the data. The performance of the laboratories was commendable, despite some variation in the reported results due to the different sequencing and bioinformatics approaches used by laboratories. The overall outcome is positive and demonstrates the critical role of the PTP in supporting the implementation and validation of SARS-CoV-2 WGS processes. The data derived from this PTP will contribute to the development of SARS-CoV-2 bioinformatic quality control (QC) and performance benchmarking for accreditation.

Genomic Epidemiology and Antimicrobial Resistance Mechanisms of Imported Typhoid in Australia.

Typhoid fever is an invasive bacterial disease of humans that disproportionately affects low- and middle-income countries. Antimicrobial resistance (AMR) has been increasingly prevalent in recent decades in Salmonella enterica serovar Typhi (S. Typhi), the causative agent of typhoid fever, limiting treatment options. In Australia, most cases of typhoid fever are imported due to travel to regions where typhoid fever is endemic. Here, all 116 isolates of S. Typhi isolated in Victoria, Australia, between 1 July 2018 and 30 June 2020, underwent whole-genome sequencing and antimicrobial susceptibility testing. Genomic data were linked to international travel data collected from routine case interviews. Travel to South Asia accounted for most cases, with 92.2% imported from seven primary countries (the top two were India, n = 87, and Pakistan, n = 12). A total of 17 S. Typhi genotypes were detected in the 2-year cohort, with 48.2% genotyped as part of global AMR lineages. Ciprofloxacin resistance was detected in two lineages, 3.3 and 4.3.1.2, all from cases with reported travel to India. Nearly all multidrug and extensively drug resistant isolates (90%) were from cases with reported travel to Pakistan in genotypes 4.3.1.1 and 4.3.1.1.P1. Extended spectrum beta-lactamases, blaCTX-M-15 and blaSHV-12, were detected in cases with travel to Pakistan and India, respectively. Linking epidemiological data with genomic studies of S. Typhi provides an opportunity to improve understanding of the emergence, spread and risk of drug-resistant S. Typhi infections and to better inform empirical treatment guidelines in returned travelers.

Evolutionary dynamics of multidrug resistant Salmonella enterica serovar 4,[5],12:i:- in Australia.

Salmonella enterica serovar 4,[5],12:i:- (Salmonella 4,[5],12:i:-) is a monophasic variant of Salmonella Typhimurium that has emerged as a global cause of multidrug resistant salmonellosis. We used Bayesian phylodynamics, genomic epidemiology, and phenotypic characterization to describe the emergence and evolution of Salmonella 4,[5],12:i:- in Australia. We show that the interruption of the genetic region surrounding the phase II flagellin, FljB, causing a monophasic phenotype, represents a stepwise evolutionary event through the accumulation of mobile resistance elements with minimal impairment to bacterial fitness. We identify three lineages with different population dynamics and discrete antimicrobial resistance profiles emerged, likely reflecting differential antimicrobial selection pressures. Two lineages are associated with travel to South-East Asia and the third lineage is endemic to Australia. Moreover antimicrobial-resistant Salmonella 4,[5],12:i- lineages efficiently infected and survived in host phagocytes and epithelial cells without eliciting significant cellular cytotoxicity, suggesting a suppression of host immune response that may facilitate the persistence of Salmonella 4,[5],12:i:-.

Proficiency testing for bacterial whole genome sequencing in assuring the quality of microbiology diagnostics in clinical and public health laboratories.

The adoption of whole genome sequencing (WGS) data over the past decade for pathogen surveillance, and decision-making for infectious diseases has rapidly transformed the landscape of clinical microbiology and public health. However, for successful transition to routine use of these techniques, it is crucial to ensure the WGS data generated meet defined quality standards for pathogen identification, typing, antimicrobial resistance detection and surveillance. Further, the ongoing development of these standards will ensure that the bioinformatic processes are capable of accurately identifying and characterising organisms of interest, and thereby facilitate the integration of WGS into routine clinical and public health laboratory setting. A pilot proficiency testing (PT) program for WGS of infectious agents was developed to facilitate widely applicable standardisation and benchmarking standards for WGS across a range of laboratories. The PT participating laboratories were required to generate WGS data from two bacterial isolates, and submit the raw data for independent bioinformatics analysis, as well as analyse the data with their own processes and answer relevant questions about the data. Overall, laboratories used a diverse range of bioinformatics tools and could generate and analyse high-quality data, either meeting or exceeding the minimum requirements. This pilot has provided valuable insight into the current state of genomics in clinical microbiology and public health laboratories across Australia. It will provide a baseline guide for the standardisation of WGS and enable the development of a PT program that allows an ongoing performance benchmark for accreditation of WGS-based test processes.

Genomics-informed responses in the elimination of COVID-19 in Victoria, Australia: an observational, genomic epidemiological study.

BACKGROUND: A cornerstone of Australia's ability to control COVID-19 has been effective border control with an extensive supervised quarantine programme. However, a rapid recrudescence of COVID-19 was observed in the state of Victoria in June, 2020. We aim to describe the genomic findings that located the source of this second wave and show the role of genomic epidemiology in the successful elimination of COVID-19 for a second time in Australia. METHODS: In this observational, genomic epidemiological study, we did genomic sequencing of all laboratory-confirmed cases of COVID-19 diagnosed in Victoria, Australia between Jan 25, 2020, and Jan 31, 2021. We did phylogenetic analyses, genomic cluster discovery, and integrated results with epidemiological data (detailed information on demographics, risk factors, and exposure) collected via interview by the Victorian Government Department of Health. Genomic transmission networks were used to group multiple genomic clusters when epidemiological and genomic data suggested they arose from a single importation event and diversified within Victoria. To identify transmission of emergent lineages between Victoria and other states or territories in Australia, all publicly available SARS-CoV-2 sequences uploaded before Feb 11, 2021, were obtained from the national sequence sharing programme AusTrakka, and epidemiological data were obtained from the submitting laboratories. We did phylodynamic analyses to estimate the growth rate, doubling time, and number of days from the first local infection to the collection of the first sequenced genome for the dominant local cluster, and compared our growth estimates to previously published estimates from a similar growth phase of lineage B.1.1.7 (also known as the Alpha variant) in the UK. FINDINGS: Between Jan 25, 2020, and Jan 31, 2021, there were 20 451 laboratory-confirmed cases of COVID-19 in Victoria, Australia, of which 15 431 were submitted for sequencing, and 11 711 met all quality control metrics and were included in our analysis. We identified 595 genomic clusters, with a median of five cases per cluster (IQR 2-11). Overall, samples from 11 503 (98·2%) of 11 711 cases clustered with another sample in Victoria, either within a genomic cluster or transmission network. Genomic analysis revealed that 10 426 cases, including 10 416 (98·4%) of 10 584 locally acquired cases, diagnosed during the second wave (between June and October, 2020) were derived from a single incursion from hotel quarantine, with the outbreak lineage (transmission network G, lineage D.2) rapidly detected in other Australian states and territories. Phylodynamic analyses indicated that the epidemic growth rate of the outbreak lineage in Victoria during the initial growth phase (samples collected between June 4 and July 9, 2020; 47·4 putative transmission events, per branch, per year [1/years; 95% credible interval 26·0-85·0]), was similar to that of other reported variants, such as B.1.1.7 in the UK (mean approximately 71·5 1/years). Strict interventions were implemented, and the outbreak lineage has not been detected in Australia since Oct 29, 2020. Subsequent cases represented independent international or interstate introductions, with limited local spread. INTERPRETATION: Our study highlights how rapid escalation of clonal outbreaks can occur from a single incursion. However, strict quarantine measures and decisive public health responses to emergent cases are effective, even with high epidemic growth rates. Real-time genomic surveillance can alter the way in which public health agencies view and respond to COVID-19 outbreaks. FUNDING: The Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.

Search and Contain: Impact of an Integrated Genomic and Epidemiological Surveillance and Response Program for Control of Carbapenemase-producing Enterobacterales.

BACKGROUND: Multiresistant organisms (MROs) pose a critical threat to public health. Population-based programs for control of MROs such as carbapenemase-producing Enterobacterales (CPE) have emerged and evaluation is needed. We assessed the feasibility and impact of a statewide CPE surveillance and response program deployed across Victoria, Australia (population 6.5 million). METHODS: A prospective multimodal intervention including active screening, carrier isolation, centralized case investigation, and comparative pathogen genomics was implemented. We analyzed trends in CPE incidence and clinical presentation, risk factors, and local transmission over the program's first 3 years (2016-2018). RESULTS: CPE case ascertainment increased over the study period to 1.42 cases/100 000 population, linked to increased screening without a concomitant rise in active clinical infections (0.45-0.60 infections/100 000 population, P = .640). KPC-2 infection decreased from 0.29 infections/100 000 population prior to intervention to 0.03 infections/100 000 population in 2018 (P = .003). Comprehensive case investigation identified instances of overseas community acquisition. Median time between isolate referral and genomic and epidemiological assessment for local transmission was 11 days (IQR, 9-14). Prospective surveillance identified numerous small transmission networks (median, 2; range, 1-19 cases), predominantly IMP and KPC, with median pairwise distance of 8 (IQR, 4-13) single nucleotide polymorphisms; low diversity between clusters of the same sequence type suggested genomic cluster definitions alone are insufficient for targeted response. CONCLUSIONS: We demonstrate the value of centralized CPE control programs to increase case ascertainment, resolve risk factors, and identify local transmission through prospective genomic and epidemiological surveillance; methodologies are transferable to low-prevalence settings and MROs globally.

Tracking the COVID-19 pandemic in Australia using genomics.

Genomic sequencing has significant potential to inform public health management for SARS-CoV-2. Here we report high-throughput genomics for SARS-CoV-2, sequencing 80% of cases in Victoria, Australia (population 6.24 million) between 6 January and 14 April 2020 (total 1,333 COVID-19 cases). We integrate epidemiological, genomic and phylodynamic data to identify clusters and impact of interventions. The global diversity of SARS-CoV-2 is represented, consistent with multiple importations. Seventy-six distinct genomic clusters were identified, including large clusters associated with social venues, healthcare and cruise ships. Sequencing sequential samples from 98 patients reveals minimal intra-patient SARS-CoV-2 genomic diversity. Phylodynamic modelling indicates a significant reduction in the effective viral reproductive number (Re) from 1.63 to 0.48 after implementing travel restrictions and physical distancing. Our data provide a concrete framework for the use of SARS-CoV-2 genomics in public health responses, including its use to rapidly identify SARS-CoV-2 transmission chains, increasingly important as social restrictions ease globally.

Mining the Methylome Reveals Extensive Diversity in Staphylococcus epidermidis Restriction Modification.

Staphylococcus epidermidis is a significant opportunistic pathogen of humans. Molecular studies in this species have been hampered by the presence of restriction-modification (RM) systems that limit introduction of foreign DNA. Here, we establish the complete genomes and methylomes for seven clinically significant, genetically diverse S. epidermidis isolates and perform the first systematic genomic analyses of the type I RM systems within both S. epidermidis and Staphylococcus aureus Our analyses revealed marked differences in the gene arrangement, chromosomal location, and movement of type I RM systems between the two species. Unlike S. aureus, S. epidermidis type I RM systems demonstrate extensive diversity even within a single genetic lineage. This is contrary to current assumptions and has important implications for approaching the genetic manipulation of S. epidermidis Using Escherichia coli plasmid artificial modification (PAM) to express S. epidermidishsdMS, we readily overcame restriction barriers in S. epidermidis and achieved electroporation efficiencies equivalent to those of modification-deficient mutants. With these functional experiments, we demonstrated how genomic data can be used to predict both the functionality of type I RM systems and the potential for a strain to be electroporation proficient. We outline an efficient approach for the genetic manipulation of S. epidermidis strains from diverse genetic backgrounds, including those that have hitherto been intractable. Additionally, we identified S. epidermidis BPH0736, a naturally restriction-defective, clinically significant, multidrug-resistant ST2 isolate, as an ideal candidate for molecular studies.IMPORTANCEStaphylococcus epidermidis is a major cause of hospital-acquired infections, especially those related to implanted medical devices. Understanding how S. epidermidis causes disease and devising ways to combat these infections have been hindered by an inability to genetically manipulate clinically significant hospital-adapted strains. Here, we provide the first comprehensive analyses of the barriers to the uptake of foreign DNA in S. epidermidis and demonstrate that these are distinct from those described for S. aureus Using these insights, we demonstrate an efficient approach for the genetic manipulation of S. epidermidis to enable the study of clinical isolates for the first time.

Genomic investigation of Staphylococcus aureus recovered from Gambian women and newborns following an oral dose of intra-partum azithromycin.

BACKGROUND: Oral azithromycin given during labour reduces carriage of bacteria responsible for neonatal sepsis, including Staphylococcus aureus. However, there is concern that this may promote drug resistance. OBJECTIVES: Here, we combine genomic and epidemiological data on S. aureus isolated from mothers and babies in a randomized intra-partum azithromycin trial (PregnAnZI) to describe bacterial population dynamics and resistance mechanisms. METHODS: Participants from both arms of the trial, who carried S. aureus in day 3 and day 28 samples post-intervention, were included. Sixty-six S. aureus isolates (from 7 mothers and 10 babies) underwent comparative genome analyses and the data were then combined with epidemiological data. Trial registration (main trial): ClinicalTrials.gov Identifier NCT01800942. RESULTS: Seven S. aureus STs were identified, with ST5 dominant (n = 40, 61.0%), followed by ST15 (n = 11, 17.0%). ST5 predominated in the placebo arm (73.0% versus 49.0%, P = 0.039) and ST15 in the azithromycin arm (27.0% versus 6.0%, P = 0.022). In azithromycin-resistant isolates, msr(A) was the main macrolide resistance gene (n = 36, 80%). Ten study participants, from both trial arms, acquired azithromycin-resistant S. aureus after initially harbouring a susceptible isolate. In nine (90%) of these cases, the acquired clone was an msr(A)-containing ST5 S. aureus. Long-read sequencing demonstrated that in ST5, msr(A) was found on an MDR plasmid. CONCLUSIONS: Our data reveal in this Gambian population the presence of a dominant clone of S. aureus harbouring plasmid-encoded azithromycin resistance, which was acquired by participants in both arms of the study. Understanding these resistance dynamics is crucial to defining the public health drug resistance impacts of azithromycin prophylaxis given during labour in Africa.

Genomics for Molecular Epidemiology and Detecting Transmission of Carbapenemase-Producing Enterobacterales in Victoria, Australia, 2012 to 2016.

Carbapenemase-producing Enterobacterales (CPE) are being increasingly reported in Australia, and integrated clinical and genomic surveillance is critical to effectively manage this threat. We sought to systematically characterize CPE in Victoria, Australia, from 2012 to 2016. Suspected CPE were referred to the state public health laboratory in Victoria, Australia, from 2012 to 2016 and examined using phenotypic, multiplex PCR and whole-genome sequencing (WGS) methods and compared with epidemiological metadata. Carbapenemase genes were detected in 361 isolates from 291 patients (30.8% of suspected CPE isolates), mostly from urine (42.1%) or screening samples (34.8%). IMP-4 (28.0% of patients), KPC-2 (25.3%), NDM (24.1%), and OXA carbapenemases (22.0%) were most common. Klebsiella pneumoniae (48.8% of patients) and Escherichia coli (26.1%) were the dominant species. Carbapenemase-inactivation method (CIM) testing reliably detected carbapenemase-positive isolates (100% sensitivity, 96.9% specificity), identifying an additional five CPE among 159 PCR-negative isolates (IMI and SME carbapenemases). When epidemiologic investigations were performed, all pairs of patients designated "highly likely" or "possible" local transmission had ≤23 pairwise single-nucleotide polymorphisms (SNPs) by genomic transmission analysis; conversely, all patient pairs designated "highly unlikely" local transmission had ≥26 pairwise SNPs. Using this proposed threshold, possible local transmission was identified involving a further 16 patients for whom epidemiologic data were unavailable. Systematic application of genomics has uncovered the emergence of polyclonal CPE as a significant threat in Australia, providing important insights to inform local public health guidelines and interventions. Using our workflow, pairwise SNP distances between CPE isolates of ≤23 SNPs suggest local transmission.

Emergence and divergence of major lineages of Shiga-toxin-producing Escherichia coli in Australia.

Shiga-toxin-producing Escherichia coli (STEC) infection is an important global cause of foodborne disease. To date however, genomics-based studies of STEC have been predominately focused upon STEC collected in the Northern Hemisphere. Here, we demonstrate the population structure of 485 STEC isolates in Australia, and show that several clonal groups (CGs) common to Australia were infrequently detected in a representative selection of contemporary STEC genomes from around the globe. Further, phylogenetic analysis demonstrated that lineage II of the global O157:H7 STEC was most prevalent in Australia, and was characterized by a frameshift mutation in flgF, resulting in the H-non-motile phenotype. Strong concordance between in silico and phenotypic serotyping was observed, along with concordance between in silico and conventional detection of stx genes. These data represent the most comprehensive STEC analysis from the Southern Hemisphere, and provide a framework for future national genomics-based surveillance of STEC in Australia.

Pleistocene-dated biogeographic barriers drove divergence within the Australo-Papuan region in a sex-specific manner: an example in a widespread Australian songbird.

Understanding how environmental change has shaped species evolution can inform predictions of how future climate change might continue to do so. Research of widespread biological systems spanning multiple climates that have been subject to environmental change can yield generalizable inferences about the neutral and adaptive processes driving lineage divergence during periods of environmental change. We contribute to the growing body of multi-locus phylogeographic studies investigating the effect of Pleistocene climate change on species evolution by focusing on a widespread Australo-Papuan songbird with several mitochondrial lineages that diverged during the Pleistocene, the grey shrike-thrush (Colluricincla harmonica). We employed multi-locus phylogenetic, population genetic and coalescent analyses to (1) assess whether nuclear genetic diversity suggests a history congruent with that based on phenotypically defined subspecies ranges, mitochondrial clade boundaries and putative biogeographical barriers, (2) estimate genetic diversity within and genetic differentiation and gene flow among regional populations and (3) estimate population divergence times. The five currently recognized subspecies of grey shrike-thrush are genetically differentiated in nuclear and mitochondrial genomes, but connected by low levels of gene flow. Divergences among these populations are concordant with recognized historical biogeographical barriers and date to the Pleistocene. Discordance in the order of population divergence events based on mitochondrial and nuclear genomes suggests a history of sex-biased gene flow and/or mitochondrial introgression at secondary contacts. This study demonstrates that climate change can impact sexes with different dispersal biology in different ways. Incongruence between population and mitochondrial trees calls for a genome-wide investigation into dispersal, mitochondrial introgression and mitonuclear evolution.

Persistent low avian malaria in a tropical species despite high community prevalence.

Malarial and other haemosporidian parasites are widespread; however, their temporal dynamics are ill-understood. Longitudinal sampling of a threatened riparian bird revealed a consistently very low prevalence over 13 years (∼5%) despite infections persisting and prevalence increasing with age. In contrast, three key species within this tropical community were highly infected (∼20-75% prevalence) and these differences were stable. Although we found novel lineages and phylogenetic structure at the local level, there was little geographic structuring within Australasia. This study suggests that malarial parasite susceptibility is determined by host factors and that species can maintain low levels despite high community prevalence.

Remodeling of pSK1 Family Plasmids and Enhanced Chlorhexidine Tolerance in a Dominant Hospital Lineage of Methicillin-Resistant Staphylococcus aureus.

Staphylococcus aureus is a significant human pathogen whose evolution and adaptation have been shaped in part by mobile genetic elements (MGEs), facilitating the global spread of extensive antimicrobial resistance. However, our understanding of the evolutionary dynamics surrounding MGEs, in particular, how changes in the structure of multidrug resistance (MDR) plasmids may influence important staphylococcal phenotypes, is incomplete. Here, we undertook a population and functional genomics study of 212 methicillin-resistant S. aureus (MRSA) sequence type 239 (ST239) isolates collected over 32 years to explore the evolution of the pSK1 family of MDR plasmids, illustrating how these plasmids have coevolved with and contributed to the successful adaptation of this persistent MRSA lineage. Using complete genomes and temporal phylogenomics, we reconstructed the evolution of the pSK1 family lineage from its emergence in the late 1970s and found that multiple structural variants have arisen. Plasmid maintenance and stability were linked to IS256- and IS257-mediated chromosomal integration and disruption of the plasmid replication machinery. Overlaying genomic comparisons with phenotypic susceptibility data for gentamicin, trimethoprim, and chlorhexidine, it appeared that pSK1 has contributed to enhanced resistance in ST239 MRSA isolates through two mechanisms: (i) acquisition of plasmid-borne resistance mechanisms increasing the rates of gentamicin resistance and reduced chlorhexidine susceptibility and (ii) changes in the plasmid configuration linked with further enhancement of chlorhexidine tolerance. While the exact mechanism of enhanced tolerance remains elusive, this research has uncovered a potential evolutionary response of ST239 MRSA to biocides, one of which may contribute to the ongoing persistence and adaptation of this lineage within health care institutions.

Co-circulation of Multidrug-resistant Shigella Among Men Who Have Sex With Men in Australia.

BACKGROUND: In urban Australia, the burden of shigellosis is either in returning travelers from shigellosis-endemic regions or in men who have sex with men (MSM). Here, we combine genomic data with comprehensive epidemiological data on sexual exposure and travel to describe the spread of multidrug-resistant Shigella lineages. METHODS: A population-level study of all cultured Shigella isolates in the state of Victoria, Australia, was undertaken from 1 January 2016 through 31 March 2018. Antimicrobial susceptibility testing, whole-genome sequencing, and bioinformatic analyses of 545 Shigella isolates were performed at the Microbiological Diagnostic Unit Public Health Laboratory. Risk factor data on travel and sexual exposure were collected through enhanced surveillance forms or by interviews. RESULTS: Rates of antimicrobial resistance were high, with 17.6% (95/541) and 50.6% (274/541) resistance to ciprofloxacin and azithromycin, respectively. There were strong associations between antimicrobial resistance, phylogeny, and epidemiology. Specifically, 2 major MSM-associated lineages were identified: a Shigellasonnei lineage (n = 159) and a Shigella flexneri 2a lineage (n = 105). Of concern, 147/159 (92.4%) of isolates within the S. sonnei MSM-associated lineage harbored mutations associated with reduced susceptibility to recommended oral antimicrobials: namely, azithromycin, trimethoprim-sulfamethoxazole, and ciprofloxacin. Long-read sequencing demonstrated global dissemination of multidrug-resistant plasmids across Shigella species and lineages, but predominantly associated with MSM isolates. CONCLUSIONS: Our contemporary data highlight the ongoing public health threat posed by resistant Shigella, both in Australia and globally. Urgent multidisciplinary public health measures are required to interrupt transmission and prevent infection.

Reconstruction of the Genomes of Drug-Resistant Pathogens for Outbreak Investigation through Metagenomic Sequencing.

Culture-independent methods that target genome fragments have shown promise in identifying certain pathogens, but the holy grail of comprehensive pathogen genome detection from microbiologically complex samples for subsequent forensic analyses remains a challenge. In the context of an investigation of a nosocomial outbreak, we used shotgun metagenomic sequencing of a human fecal sample and a neural network algorithm based on tetranucleotide frequency profiling to reconstruct microbial genomes and tested the same approach using rectal swabs from a second patient. The approach rapidly and readily detected the genome of Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae in the patient fecal specimen and in the rectal swab sample, achieving a level of strain resolution that was sufficient for confident transmission inference during a highly clonal outbreak. The analysis also detected previously unrecognized colonization of the patient by vancomycin-resistant Enterococcus faecium, another multidrug-resistant bacterium.IMPORTANCE The study results reported here perfectly demonstrate the power and promise of clinical metagenomics to recover genome sequences of important drug-resistant bacteria and to rapidly provide rich data that inform outbreak investigations and treatment decisions, independently of the need to culture the organisms.

Comprehensive antibiotic-linked mutation assessment by resistance mutation sequencing (RM-seq).

Mutation acquisition is a major mechanism of bacterial antibiotic resistance that remains insufficiently characterised. Here we present RM-seq, a new amplicon-based deep sequencing workflow based on a molecular barcoding technique adapted from Low Error Amplicon sequencing (LEA-seq). RM-seq allows detection and functional assessment of mutational resistance at high throughput from mixed bacterial populations. The sensitive detection of very low-frequency resistant sub-populations permits characterisation of antibiotic-linked mutational repertoires in vitro and detection of rare resistant populations during infections. Accurate quantification of resistance mutations enables phenotypic screening of mutations conferring pleiotropic phenotypes such as in vivo persistence, collateral sensitivity or cross-resistance. RM-seq will facilitate comprehensive detection, characterisation and surveillance of resistant bacterial populations ( https://github.com/rguerillot/RM-seq ).

Global spread of three multidrug-resistant lineages of Staphylococcus epidermidis.

Staphylococcus epidermidis is a conspicuous member of the human microbiome, widely present on healthy skin. Here we show that S. epidermidis has also evolved to become a formidable nosocomial pathogen. Using genomics, we reveal that three multidrug-resistant, hospital-adapted lineages of S. epidermidis (two ST2 and one ST23) have emerged in recent decades and spread globally. These lineages are resistant to rifampicin through acquisition of specific rpoB mutations that have become fixed in the populations. Analysis of isolates from 96 institutions in 24 countries identified dual D471E and I527M RpoB substitutions to be the most common cause of rifampicin resistance in S. epidermidis, accounting for 86.6% of mutations. Furthermore, we reveal that the D471E and I527M combination occurs almost exclusively in isolates from the ST2 and ST23 lineages. By breaching lineage-specific DNA methylation restriction modification barriers and then performing site-specific mutagenesis, we show that these rpoB mutations not only confer rifampicin resistance, but also reduce susceptibility to the last-line glycopeptide antibiotics, vancomycin and teicoplanin. Our study has uncovered the previously unrecognized international spread of a near pan-drug-resistant opportunistic pathogen, identifiable by a rifampicin-resistant phenotype. It is possible that hospital practices, such as antibiotic monotherapy utilizing rifampicin-impregnated medical devices, have driven the evolution of this organism, once trivialized as a contaminant, towards potentially incurable infections.